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We consider the acoustic propagation in lined flow duct with a purely reactive 
impedance at the wall. This reacting liner has the capability to reduce the speed 
of sound, and thus to enhance the interaction between the acoustic propagation and 
the low Mach number flow (M ~ 0.3). At the lower frequencies, there are typically 
4 acoustic or hydrodynamic propagating modes, with 3 of them propagating in the 
direction of the flow. Above a critical frequency, there are only 2 propagating modes 
that all propagate in the direction of the flow. From the exact 2D formulation an 
approximate ID model is developed to study the scattering of acoustic waves in a 
straight duct with varying wall impedance. This simple system, with a uniform flow 
and with a non-uniform liner impedance at the wall, permits to study the scattering 
between regions with different waves characteristics. Several situations are charac¬ 
terized to show the importance of negative energy waves, strong interactions between 
acoustic and hydrodynamic modes or asymmetric scattering. 
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I. INTRODUCTION 


Acoustic liners are widely used to reduce the sound transmission in ducts with flow 
with applications in household appliances, ventilation systems in vehicles and buildings, IC- 
engines, power plants, aircraft engine. The mitigation due to these liners is based on two 
principles that are generally mixed. The first action of the liners is to dissipate acoustical 
energy by visco-thermal losses or by exchange of energy between the acoustical held and 
flow, like in the vicinity of a hole in a perforated plate with grazing flow. The second type 
of action is the scattering of the acoustical waves by the changes of acoustical impedance 
occurring for instance at the entrance and at the exit of the liners. This paper focusses on the 
second type of action called reacting effects and disregards the first type called dissipative 
effects. To do this, a waveguide with an acoustically treated wall is studied and the wall is 
considered as locally reacting and without dissipation. When the liner consists of cavities 
mounted flush to the wall (like small closed tubes in the present case), those cavities act as 
springs in the low frequencies limit. Then, the speed of sound is determined by the square 
root of the ratio between the isentropic bulk modulus (which is a measure for the stiffness 
of the fluid) and the mass density. The presence of small cavities decreases the effective 
stiffness and, consequently, the speed of sound. The propagative acoustical waves in such 
systems are called ’’slow sound”. Recently, slow acoustic waves have attracted attention for 
the potential to design new acoustic devices such as metamaterials. They have been studied 
both in sonic crystal^ and in one dimensional systeirpEl. The originality of the present study 
is to introduce a mean flow with a velocity of the same order as the effective speed of sound. 
When the flow velocity is smaller than the speed of sound the regime is called ’’subsonic” 
and ’’supersonic” on the other case. In the subsonic regime, it will be shown that 4 modes 
propagate at low frequencies (wavelength much smaller than the transverse dimensions of the 
waveguide). Two of these modes correspond to classical acoustical waves in both direction. 
The two other modes do not exist without flow and are thus called HydroDynamic (HD) 
modes in the following. One of these HD modes has a group velocity and a phase velocity 
in opposite direction. The second HD mode is a Negative Energy Wave (NEW). Globally, 
among the 4 modes that propagate in the subsonic regime, 3 of the modes propagate in 
the flow direction while one of the acoustical modes propagates against the flow. In the 
supersonic regime, only 2 waves can propagate and they are in the flow direction. The 
problem that we consider is close to the response of fluid loaded finite plates with mean 
flov3® but it leads to a simpler analysis of interesting behaviors. 


The plan of the paper is as follows. The section |TT] of this paper is devoted to the 
characterization of the modes propagating in the low frequencies limit in a 2D duct. It 
will be shown that an energy flux conservation can be written in this case. The section 


HI describes a approximate ID model where the effects in the transverse direction of the 


duct are taken into account by averaging. Albeit very simple, this ID model has the same 
richness of behavior as the 2D model. In particular, as in the 2D model, an energy flux 
conservation is obtained and a NEW is present. In the section IV, the ID model is used 


to calculate the scattering proprieties of an increase or a decrease in the wall impedance. 
The transonic cases (from supersonic to subsonic and vice versa) are of particular interest 
because of the conversion of acoustical waves into HD modes. A local transonic increase of 
the compliance is also studied in the section |V] and shows an interesting propriety of total 
transmission in flow direction and of no transmission in the opposite direction corresponding 
to an ’’acoustical diode”. 
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II. SOUND PROPAGATION IN A 2D DUCT WITH FLOW AND 
COMPLIANT WALL 



We consider the sound propagation in a 2D channel with a uniform flow, see Fig. 1. The 
lower wall is rigid. The upper wall is compliant and composed of small tubes of variable 
lengths. All parameters are nondimensionalized in the standard way to simplify the notation. 
Velocities are nondimensionalized by the speed of sound Cq, so that the uniform mean velocity 
becomes the steady flow Mach number M. Distances are nondimensionalized by the height 
of the channel H, time by H/cq, and pressure by PqCq where po is the mean density. The 
dimensionless equations governing the acoustic motion are then 

Dtp=-V.\ (1) 

A V = -Vp , (2) 

where p is the pressure, v is the velocity and Dt = dt + Md^ is the convective derivative. 
Next, the equations are written in term of the acoustic velocity potential (v = Vp). Eq. 
Q leads to p = —Dt (f and Eq. Q leads to the classical convected wave equation; 

AV - vV = 0. (3) 

On the rigid wall (p = 0), the boundary condition is dyp = 0. On the compliant wall, we 
use the so-called ’Tngard-Myers Condition”®. This condition states that the pressure and 
the transverse displacement rj {Dtp = v = dyp) are continuous at the wall which leads to 
V = Dt{C{x)p) where C{x) is the compliance of the wall (ratio of transverse displacement 
over the pressure). Hence the boundary condition at the wall p = 1 is written: 

dy(f = -Dt {C{x)Dt(p). (4) 

The compliance of the closed tubes of length b{x) at p = 1 is given by C{x) = 
a tein{b{x)u)/u where a is the percentage of open area (POA, ratio between the surface 
of the tubes and the total surface). In the very low frequencies limit {ub -C 1), the com- 
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pliance is simply equal to the length of the tubes b{x) multiply by the POA. It means that 
in this limit, the closed tubes act like springs of stiffness 1/ab. To simplify the notation in 
the following a is suppose to be equal to unity (it could be integrated very easily if it differs 
signihcantly from unity) and the problem can be written globally as 

( D‘fip - VV = 0 

< dyip = 0 at ?/ = 0 (5) 

ydyp) =-Dt{b{x)Dtip) BXy = l. 

The impedance boundary condition with uniform flow is questionable and more advanced 
models exislP. The Ingard-Myers condition has been used here for simplicity. In the low 
frequencies limit, more complex models have been tested without qualitative changes in the 
behavior^. Furthermore, it could be noted that the low frequencies limit used here, leads 
to a ’’well posed” problem in the sense given by Bramblej^. 

A. Dispersion equation in the 2D problem 

For uniform compliance b, the solution is searched under the form ip = A cosh (a?/) exp(i(—a;f+ 
kx)) where = k'^ — and = a; — Mk. This leads to the dispersion equation: 

V{uj, k) = a tanh(Q;)-^^ = 0 (6) 

UJ 

which, in the very frequency limit, becomes: 

'D{u:, k) = a tanh(a) — 6 = 0. (7) 

Without flow in the very low frequency limit, the dispersion equation ([^ can be simplihed 
to = (1 + b)u‘^. The phase velocity 


u 1 

k a/1 + 6 


( 8 ) 


is always smaller than 1, meaning that Cb is smaller than the speed of sound in free space. 
Thus the acoustic wave propagation can be signihcantly slowed down in a duct with a wall 
which reacts locally like a spring. The phase velocity of this ’’slow sound” can be decreased 
to become of the order of the how velocity in the duct. In this case, dramatic ehects of a 
how with moderate Mach number (M ~ 0.3) are expected. 

The dispersion curves are displayed on Fig. 2 and show the ehect of the how. For 
w < iXmax, Eq. (g has 4 real solutions corresponding to propagating modes. Two roots, 
labeled S and A~ on Fig. 2, approach to each other when uj —>■ oJmax- They coalesce for a 
frequency cOmax above which they no longer exist as real roots i.e. as propagating waves. It 
is shown in Appendix]^ that when M > Cb there are only two propagating modes whatever 
the frequency (i.e. Umax = 0). 

In the subsonic regime, i.e. when M < Cb and u < Umax, 2 of the 4 modes have a 
vanishing wavenumber when a; —)■ 0. These solutions are called acoustic and, in the low 
frequencies limit, they propagate in both directions with the speed of sound Cb corrected 
by the convective effects. The 2 other solutions do not exist without flow and are called 
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Figure 2: (Color online) Propagative wave numbers k as a 
function of the frequency oj for 6 = 4 and M = 0.3. The 2 
embedded boxes give the direction of the waves 
propagation. 


HydroDynamic (HD) modes. The solution called S in Fig. 2 has a negative phase velocity 
c# = u/k but a positive group velocity Cg = du/dk and thus propagates in the flow direction. 
The last solution called N has both positive phase and group velocities and propagates in 
the flow direction. It can be seen from Fig. 12 in Appendix]^ that D < 0 for this wave and it 
will be seen below that it corresponds to a Negative Energy Wave (NEW). In the supersonic 
regime, i.e. when M > ct or u > Umax, only the A+ and the N waves can propagate. 


In summary, in the subsonic case 4 waves propagate. Three of them propagate in the 
flow direction (A+, S, N) and one propagates against the flow (A“). In the supersonic case, 
only 2 waves (A+, N) can propagate and they are in the flow direction. A NEW is always 
present. If M > 1/a/1 + b, we are always in the supersonic case. If M < 1/a/1 + b, the 
transition from subsonic to supersonic can been reach either by increasing a; at a given b or 
by increasing 6 at a given u. This last possibility will be used in the Section on scattering 
(Sec. IV and 0 . 


B. Energy flux conservation of slow sound waves with flow 


Thereafter the problem is studied in the frequency domain (convention e“‘“*) where dt = 
—lu and = — ica + Mdx- In order to And an ’’energy like” equation, the Eq. (|^ is 
classically multiplied by ip (the complex conjugate of p) and is integrated on the cross 
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section to yield: 


( / ifiidl 




- Dl^) ) = 


dx - MtpD^cp) d|/^ + [{^dyLpf^ = 

[ {ipdx(p - MlpD^if) dy 


-MLp{x,l)b{x)D^^p{x,l)) = 0, 


(9) 


where the relation 

Sm {gD^{f{x)D^g)) = (^>ni {Mgf\x)D^g)) (10) 

valid for any function g and any real function / had been used. Thus the quantity 


J = 9=ni 


g:>{d^^ - MD^ip) dy 


-Mp}{x, l)b{x)D^p}{x, 1) 


( 11 ) 


is conserved along x. This expression is identical to the expression of the energy flux proposed 
by Mohrin^l^. This energy flux can be computed for each mode m (normalized by its value 
at y = 1): ipm = cosh.{amy)G^^'^^/ cosh(Q;m) where am is one of the solutions of the dispersion 
equation, and flm = oj — Mkm '■ 


T — i (u I /\/rn + 20;^ , ]\jun 

Jm ( (^m “1“ 1 2/ \ Mhi.ly 

4a^cosh [a,n) 


The value of Jm is displayed on Fig. 3 for the 4 propagating modes. The modes and 
S', that propagate in the flow direction, have positive energy fluxes. The modes A~, that 
propagates against the flow, has negative energy flux. The mode iV, that propagates in the 
flow direction, has a negative energy flux. It means that this mode is a NEW®. This last 
property will have important consequences in the results presented afterwards. 


III. ID MODEL 
A. ID approximation 

In order to simplify the analysis of the problem, we are looking for a ID model that 
conserved the main proprieties of the 2D problem: The dispersion relation has to give the 
same number of propagating modes as the 2D model and a conserved energy flux has to be 
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the frequency u for 6 = 4 and M = 0.3. 


defined. For that, we integrate the 2D equation ([^ along y and we get the exact expression: 

Df ^ d?/^ (p d?/^ - dy^{x, 1) = 0 (12) 

which is associated to the boundary conditions in ([^. A simplihcation can be achieved if 
we now assume that the y derivative of (p at the compliant wall can be written: 


dyip{x, 1) = aiV {x) + a 2 F{x) 


(13) 


where Oi and 02 are two real constants and where V and F are dehned as the two functions 
appearing in (12) and in the boundary condition at j/ = 1 in Eq. (|^: 


V(x)=(p(x,l) and F(x) 



(pdy. 


This leads to the system of two coupled ODEs 


D^F - dlF = aiV + asF 
DtihDtV) = -(aiE + aaF) 


(14) 

(15) 


The real constants ai and 02 can be chosen freely. For instance, for a parabolic approxi¬ 
mation such as ip = Cl + C 2 y‘^, the constants are Oi = —02 = 3. 
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B. Dispersion relation 


When b is constant, looking for a solution under the form exp(i(—cut + kx)), leads to the 
dispersion equation expressed as a function of the frequency in the moving frame hi = u:—Mk\ 


ip? — k'^ + 02 ) — oi) + 0102 = 0 (16) 

The solutions in term of u versus k are plot on Fig. 4 and compared to the results of the 
2D model. The agreement between the two model is good when the coefficient oi = —02 is 
chosen in such a way that the value of k when a; —?■ 0 of the modes N and S are closed in 
ID and 2D model. When oi = — 02 , the wavenumbers for a; —?• 0 are: 




+ b u! 

1 ± M \/l + b 


and /civ 5 = 


/oi(l-M2(l + 6)) 

6M2(1 - M2) 


(17) 



Figure 4: (Color online) Propagative wave numbers k as a 
function of the frequency io for oi = —02 = 4, 6 = 4 and 
M = 0.3. The solid lines represent the solutions of the ID 
model while the dashed lines represent the solution of the 
2D model (see Fig. 2). 


C. Evolution equations 

A set of first order evolution equations can be derived from the Eqs. and ^ by 
introducing G and W with G = (a; + i(l — M)dx)F and (m + iMd.j;)W = —(aiV + 02 ^). In 











vectorial notation, the evolution equation is: 


and 


-i5,,X 


QX where 


/F\ 

G 

V 

\wj 


-u/il-M) 1/(1-M) 0 0 

02 /( 1 + M) a;/(l + M) ai/(l + M) 0 
0 0 uj/M 1/(M6) 

a 2 /M 0 ai/M uj/M 


(18) 


The eigenvalues of the matrix Q are the four k^, solutions of the dispersion equation 
(16) and the eigenvectors X^ give a relation for each modes between the mean value of 
the velocity potential over the section and its value at the wall Vm- Note that, at low 
frequencies, the A~^ and A~ modes are quasi plane while the S and N are more localized 
along the compliant wall. 


D. Energy flux conservation in the ID model 


To obtain an energy flux conservation, we multiply the Eq. (14) by F and Eq. (15) by 
V and we make use of the relation (10) : 


Thus the quantity 


(Om (Fd^F - MFD^F)) = (FV) 

a, (3m {MbVD^V)) = -023111 (VF) 


I = 3m G(a^F - MD„F) + (/m hVD^V 


(19) 

( 20 ) 

( 21 ) 


is conserved along x. It can be noticed that if Ui = — 02 , the ID energy flux conservation 
(21) becomes very similar to the exact energy flux conservation in 2D, see equation In 
this case, the energy flux of any mode (m = y4+, A~, S', X) is given by: 


Im = (fcni(l - M^) + Mu:)\Fm? + M6(a; - krnM)\V^\\ (22) 

As in the 2D case, the modes A"*" and S have positive energy fluxes while the modes A~ and 
N have negative energy flux. The mode X, propagating to the right, is thus a NEW, as in 
the 2D model. 

The ID model reproduces correctly all the main physical ingredients (dispersion and 
energy flux conservation) that are present in the 2D model. This model will be used in the 
next section to study the scattering induced by a change in the wall compliance in the next 
section. 
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Figure 5: (Color online) Scattering problems by a change 
in the wall compliance. The propagative waves are given 
for the 4 cases considered. Case 1: subsonic everywhere 
(The scattering matrix S is 4 x 4). Case 2: transcriticai 
variation of the compliance from supersonic to subsonic (S 
is 3 X 3j. Case 2: transcriticai variation of the compliance 
from subsonic to supersonic (S is 3 x 3). Case 4: 
supersonic everywhere fS is 2 x 2). 


IV. SCATTERING BY A CHANGE IN THE WALL COMPLIANCE 

We consider the problem defined in Fig. 5: The compliance of the wall is changing 
around x = 0 from the value bi (x < 0 ) to a value 62 > 0 ), the flow being in the positive 

x-direction. The 4 cases indicated on Fig. 5 will be considered. 


A. Subsonic case (case 1) 

In the case 1, the problem is subsonic everywhere. Upstream, at left, there are 3 incoming 
and 1 outgoing waves. Downstream, at right, there are 1 incoming and 3 outgoing waves. 
The upstream and downstream propagative field can be described by: 

X,-(x) = aj±f 

+s,Xfe‘'^^ + ajXf e'^^r- (23) 

where j = 1 or 2 labels the region, the hat indicates that the modes have been normalized 
such that their energy flux is 1 for the modes A+ and S and —1 for the modes A~ and N 

(X- = x-/^7j). 

The effect of the compliance variation is described by the scattering matrix linking the 4 
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outgoing waves B to the 4 incoming waves A: 



/4\ 


/a+\ 

B = SA where B = 

n2 

and A = 

ni 


S2 


Si 






(24) 


where the coefficients of S are classically given by Smn with n the incident mode and m the 
outgoing mode. 

The classical unitary relation!^ for a conservative system is replaced in our case (see 
Appendix [B|) by 


S^JS = J 


(25) 


where J = diag(l, —1, 1, 1) and the superscript T denotes the transpose operation. 

Even if the scattering matrix can be easily computed for a continuous variation of h{x) by 
numerical integration, for the sake of simplicity only the results for discontinuous variation 
of 6 at a; = 0 will be presented. It can be seen from Eq. (18) that the functions F, G, V and 


W are continuous when b is discontinuous while the slope of V is discontinuous at a; = 0. 
The continuity of X at x = 0 can be written in vectorial form, separating the incoming and 
the outgoing waves: 


^A+ 

^2 5 ^2 5 ^2 5 


Vo 


1 5 5 5 ^2 

' ’v' 

V, 


The scattering matrix is then computed by: 


/at\ 

n2 

S2 

V«r/ 

Ui 

Si 

w 


(26) 


S = Vo' V, 


(27) 


As an example of the scattering matrix elements, the value of the outgoing waves when the 
wave Af is incident are plotted in Fig. 6. 

It can be seen that the wave Af is mainly transmitted on A^ and some acoustical re¬ 
flection on A^ occurs. The acoustical transmission and reflection are nearly constant up to 
the value Umax = 0.1228 where the propagation becomes sonic in the tube with the larger 
b. We can dehne the acoustical transmission coefficients by = S^+^+ and the acoustical 
reflection coefficients by R\ = 5'^-^+. It can be seen from Fig. 6 that an energy-like conser¬ 
vation for the acoustical waves can be written: |T^p-|- is close and always smaller than 
1. There is some conversion from the acoustical modes to the HD modes N and S. This 
conversion increases linearly from 0 at cu = 0 but those modes are such that their energies 
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Figure 6 : (Color online) Value of the outgoing waves 
when the wave Af is incident (a)( = 1 j for M = 0.3, 

6 i = 4, 62 = 1 and ai = —02 = 4. The 4 curves represent 

the 4 outgoing waves: —; +++: N 2 , 000 : S 2 and - : 

. The thin dashed line represent the value of 
I T^ P + I RX I ^ • The embedded figure is a zoom around 1. 


are opposite to fulfill the exact energy conservation, from Eq. (21): 

\TX\^ + \RX\"+\S..,at\"-\Sn 


n2a-[ 


= 1. 


When HD modes Ni or Si are incident (not displayed), they are mainly transmitted on the 
same mode N 2 resp. S 2 and some extra transmission occurs on S 2 resp. N 2 . The conversion 
to acoustical modes is small. 

The overall picture of the subsonic case is that both the acoustical modes and the HD 
modes are rather independent. Some small conversions exist between those two families of 
modes. When the problem is near transonic, the coupling between the different kind of 
modes becomes larger. 


B. Transonic case (case 2 ) 

In the transonic case 2, on the upstream side M > 1/a/1 + bi and the propagation is 
supersonic whatever u. In this case, the modes Si and are no longer propagative but 
they are transformed into 2 evanescent modes that are complex conjugate: E~^ and E~. The 
E^ mode is dehned such as it decreases when x increases {‘Am{kE+) >0). In this transonic 
case, 2 incoming waves are present upstream and 1 incoming and 3 outgoing waves are 
present downstream. Therefore the scattering matrix S is now a 3 x 3 matrix. To apply 
the continuity of X at x = 0, it is necessary to take into account the evanescent mode 
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that decays in the —x direction {E ). The output matrix Vq is transformed into Vq = 
X^, Xf, —Xf ] while the input matrix V| is reduced to V| = [X^"*^, X]^, —X^ ]. The 
scattering matrix is obtained from Vq^ V| by removing the last line linked to the evanescent 
mode. The coefficient of the scattering matrix are now complex numbers and the absolute 
values of 3 of these coefficients (when A 2 is incident) are shown in Fig. 7. 



waves when the wave A 2 is incident (a^ = Ij for M = 0.3, 
61 = 12, 62 = 4 and ai = —02 = 4. The 3 curves are linked 
to the 3 outgoing waves: —; +++: N 2 and 000 : 82 - 


The hrst striking point on the Fig. 7, is the divergence of 2 of the curves at a; —)■ 0. When 
a; —)■ 0, the energy flux, see Eq. (22), for all the acoustical modes go linearly to 0 while the 
energy of HD modes (with a non-zero value of fc at a; = 0) do not go to zero in the subsonic 
region. In the supersonic region, the energy of the N mode goes to 0 when ca —)■ 0. To 
ensure continuity in X, the amplitude of some of the mode coefficients has to go to inhnity 
(like while the amplitudes of the eigenvectors go to 0 due to the normalization. 

When an acoustical wave is incident from the upstream side, its transmission is closed to 
1 (not displayed). Nevertheless, two HD modes with opposite energy flux are created. When 
the mode N is incident upstream, this wave is mainly transmitted with an amplitude larger 
than 1 due to the negative energy characteristic of the wave. The S wave is also created but 
the conversion into acoustical wave is weak. Interestingly, when an acoustical wave is send 
downstream H 7 , see Fig. 7, it is mainly converted into S and N waves and the reflection 
on the acoustical wave is weak (the absolute value of the acoustical reflection coefficient 
is of the order of 0.15). Thus, most of the incident acoustical energy had been transferred 
to the HD modes. This fact is also illustrated in Fig. 8 where a temporal simulation of the 
Eqs. (14) is given. It can be also remarked in this hgure that the group velocity of the 2 
HD modes are close (they are equal when ca —)■ 0) and much smaller than the group velocity 
of A'f (resp. 0.170, 0.185 and 0.747 in the present case). 
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Figure 8 : (Color online) Temporal simulation of the 
variable F in the case 2. The different colors represent the 
absolute value of F in logarithmic scale for M = 0.3, 
bi = 12 , 62 = 4 and ai = —02 = 4. An A 2 pulse is send at 
t = 0 with a central frequency Us = 0.015. The slope of the 
black dashed lines is determined by calculating the group 
velocity of each modes. 


The region x < 0 is a region which cannot be excited from the outside, although waves 
can escape from it. In particular, NEW can escape from this region. In this sense, it can be 
seen as an acoustical analogous of a white hole in general relativitj^. 


C. Transonic case (case 3 ) 

The transonic case 3, can be treated with a method similar to the case 2 except that 
the evanescent wave that had to be taken into account is Ei^. The output matrix Vq is 
transformed into Vq = [X^^,X^,Xf'^, —] while the input matrix V| is reduced to 
V| = [X(^^,X(^, —Xf]. The scattering matrix is obtained from Vq^ V| by removing the 3th 
line linked to the evanescent mode. 

When an acoustical mode is incident upstream, see Fig. 9, its transmission is again 
close to 1 with a small acoustical reflection. As a N wave is created, the acoustical energy 
increases and is close and always greater than 1. When an HD modes N or S 

is incident (not displayed), it creates a transmitted N mode and a large part is reflected as 
an acoustical mode. Only a small part is transmitted as an acoustical mode. 

The region x > 0 is a region from which no wave can escape. This can be seen as a dumb 
hole, i.e. an acoustic analogue of a black hold^. The presence in this new analogue system 
of effective horizons opens up new possibilities to explore the black hole evaporation with 
experiment^. 
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Figure 9: (Color online) Absolute value of the outgoing 
waves when the wave Af is incident (af = Ij for M = 0.3, 
5i = 4, 62 = 12 and ai = — 02 = 4. The 3 curves are linked 

to the 3 outgoing waves. —; A 2 , +++'■ N 2 and - ; A'^. 

The thin dashed line represent the value of 

The embedded figure is a zoom around 1. 


D. Supersonic case (case 4 ) 

When the problem is supersonic everywhere 2 incoming waves are present upstream and 
2 outgoing waves are present downstream (Fig. 5). Then there is no reflection and the 
scattering matrix is reduced to a 2 x 2 matrix that can be computed in the same way as 
previously, taking into account an evanescent mode on both sides of the discontinuity. The 
waves are mainly transmitted. Due to the characteristic of the negative energy wave, the 
transmission of the waves is always larger than 1 . 


V. SCATTERING BY A LOCAL INCREASE IN THE WALL 
COMPLIANCE 

A local increase in the wall compliance is depicted in Fig. 10. This configuration is 
computed as previously: the continuity of X is applied at x = 0 and at x = L, the prop¬ 
agation of the 2 modes and N 2 is taken into account between x = 0 and x = L and 
2 decreasing evanescent modes are present on each side of the compliance bump (E^ and 
E 2 ). There are 8 unknowns C = [og , ng, Sg, a^, at, ^ 2 , et, , ]'^ and four input values: 
A = [a^, rii, Si Og ]'^. From the 8 continuity relations, a vectorial relation can be written as 
McC = MaA. The global 4x4 S-matrix is composed of the first four lines of the matrix 
The absolute values of coefficients Tt = S++, R\ = S-+,T7 = S-- and 
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Ra = Sa+a- are displayed in Fig. 11. 
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Figure 11: (Color online) Absolute value of acoustical 
transmission and reflection coefficients in the Bow direction 
and R\ and against the flow and R~^ for M = 0.3, 
bi = 12, 62 = 4, ai = —02 = 4 and L = 5. 


The acoustic transmission in the flow direction is close to 1 while acoustic transmission 
against the flow is close to 0 because no wave can propagate against the flow. In this low 
frequency range, some acoustic is transmitted by the evanescent modes (tunneling effect). 
This system has been completely asymmetrised by the flow and act as an ” acoustical diode” 
for a large range of frequencies. 


VI. CONCLUDING REMARKS 

We have shown that the propagation of slow sound with flow at moderate Mach number 
have interesting and new properties. With such a system, it is possible to have subsonic 
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and supersonic propagation and to make transition (soft or abrupt) from one regime to the 
other. The scattering properties of those transitions are very similar to what happens to 
light near a white or a black hole. This analogy is useful in two ways. For instance, an 
acoustical analogue of a ’’black hole laser”^can be studied as an inverse of the work done in 
Section 4 (Supersonic —)■ Subsonic —)■ Supersonic). On the other hand, this new acoustical 
analogy open opportunities to do simple experiments on these subjects. 
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Appendix A: Dispersion relation in G vs fc 

The dispersion curves are displayed on Fig. 12 in term of G vs k. The thick continuous 
curves (G-curves) represent the Eq. (|^; \/k‘^ — tanh(-\/fc2 — — tan( 6 a;)G^/a; = 0. At 

low frequencies, this curve depends only on b. The straight line represents G = cu — kM 
and the solutions are found at the intersection of the G-curves and of the straight line. In 
the displayed case, Eq. ([^ has 4 real solutions. When ui increases at hxed M (parallel 
translation of the straight line toward larger G), two roots, labeled S and A~ on Fig. 12, 
become closer and closer. They coalesce for the frequency Umax- The same phenomenon 
occurs when M increases at fixed u (rotation of the straight line toward larger negative 
slopes). 

It can be also seen in Fig. 12 that if the slope of the straight line (given by M) is 
larger than the slope of the dispersion relation at the origin (dG/dfc = Cb where Cb is the 
’’slow sound” velocity given by Eq. (|^), represented by dashed lines in the hgure, only 
two solutions can exist whatever u. In the ID model, the slope at the origin is equal to 
±1 /a/1 — a 2 b/ai. In order to have the same slope at the origin in the ID and 2D model, we 
must have 02 = —oi. 

The Q vs k representation is also interesting because it allows the determination of the 
group velocity in the moving frame linked to the flow = dD/dfc. In our case, 3 solutions 
have a negative (A“, S, N) and one has a positive (^^)- Depending on the sign of the 
curvature of the D-curve, the group velocity at high k can be larger or smaller than the group 
velocity at low k. These cases are usually referred to in the literature as ’’superluminal” and 
’’subluminal” dispersion relation^. The slow sound analogy has a ’’subluminal” dispersion 
relation. 


Appendix B: Modified unitary relation 


When the wave is taken under the form (23), the energy flux conservation (21) between 
the regions 1 and 2 can be written: 


~ W‘2\^\ + k2n ~ (Bi) 
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Figure 12: (Color online) Solutions of the dispersion 
equation Eq. ([^ in term off! versus k for 6 = 4. The 
solution in term of of u versus k can be found at the 
intersections between these fl-curves and the straight line 
Q = Lo — kM with uj = 0.05 and M = 0.3. The thin dashed 
lines represent the slope at the origin and the thick dashed 


lines is the fl-curve of the ID model (Eq. (16) with 
0.2 = -oi = 4j. 


Splitting the incoming and the outgoing modes leads to: 

i«^r-imri + kiri + i4r = 

- l^2n + |s2n + (B2) 


which can be written: 


— T — T 

A JA = B JB 


(B3) 

where the vectors A and B are given in Eq. (24) and J = diag(l, —1, 1, 1). Using the 


dehnition of the scattering matrix, it can be written: 


— 'T — 7 ^— T' 

A JA = A S JSA 


(B4) 


This relation, valid whatever A, leads to the Eq. (25) 
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